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Abstract 

The design of a linear parameter varying (LPV) con- 
troller for an aircraft at actuator failure cases is pre- 
sented. The controller synthesis for actuator fail- 
ure cases is formulated into linear matrix inequality 
(LMI) optimizations based on an estimated failure 
parameter with pre-defined estimation error bounds. 
The inherent conservatism of an LPV control syn- 
thesis methodology is reduced using a scaling factor 
on the uncertainty block which represents estimated 
parameter uncertainties. The fault parameter is esti- 
mated using the two-stage Kalman filter. The sim- 
ulation results of the designed LPV controller for 
a HiMAT (Highly Maneuverable Aircraft Technol- 
ogy) vehicle with the on-line estimator show that 
the desired performance and robustness objectives 
are achieved for actuator failure cases. 


1 Introduction 

A fault tolerant control (FTC) system is one of the 
important research areas for a single aircraft con- 
trol upset prevention. A FTC system requires to 
detect sensor/actuator faults and to reconfigure con- 
trol laws according to faults. Also, fault detection 
and isolation (FDI) mechanism in a FTC system 
should properly integrate with a fault tolerant con- 
troller which can be designed as a reconfigurable 
controller. A fault tolerance controller may degrade 
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performance but should prevent fetal failure of a a\ i- 
ation system for fault accommodation. 

Numerous research in the field of a FTC system 
has been done. 1 * ^ 3 In Ref.[ 2,3 ], a fault tolerant con- 
troller is designed based on a priori information on 
possible faults which can be modeled as a function 
of parameters. When possible faults are modeled as 
a linear parameter varying (LPV) system, LPV con- 
trol synthesis is naturally fit into designing a fault 
tolerant controller. An LPV control synthesis is for- 
mulated into a linear matrix inequality (LMI) opti- 
mization based on an LPV system whose state-space 
matrices are functions of a scheduling parameter vec- 
tor. A scheduling parameter vector should be mea- 
surable at current time and its rate is bounded. 4 5 
The LPV controller robustly stabilize the close-loop 
system and achieve the desired performance level 
over the entire parameter space. 4,5 The LPV syn- 
thesis methodology has been successfully applied to 
synthesis controllers for the pitch-axis missile au- 
topilots, 6,7 F-14 aircraft lateral-directional axis dur- 
ing powered approach, 8 turbofan engines 9 and F-16 
aircraft. 10 

In this paper, a fault tolerant controller is de- 
signed using an LPV control synthesis methodology. 
However, the fault parameter can not be directly 
measured via sensors at current time, since the fault 
parameter is not a physical quantity. Using a FDI 
mechanism, the fault parameters can be estimated 
with a bounded error. 11 Estimation errors always 
exist in a real control system due to unmodeled dy- 
namics, disturbances, and measurement sensor er- 
rors. Thus, parameter estimation errors should be 
considered in a fault-tolerant control design process. 

In this paper, fault estimation errors and model- 
ing uncertainties are represented by an uncertainty 
block in the construct of a LPV controller. The 
structure of an uncertainty block is included in a 
robust LPV control synthesis methodology. The ro- 
bust LPV control synthesis is formulated into two 
LMI optimization problems which can be solved by 


1 

American Institute of Aeronautics and Astronautics 



an iteration method. 4 

In this paper, actuator failures are modeled as 
an LPV system as functions of actuator effective- 
ness parameters. 11 These parameters are estimated 
as biases using an augmented Kalman filter. A set 
of covariance-dependent forgetting factors is intro- 
duced into the filtering algorithm. As a result, the 
change in the actuator effectiveness is accentuated to 
help achieve a more accurate estimate more rapidly. 
The Hoo bounds on parameter estimation errors are 
assessed through simulations, which are then used 
as bounds of real parameter uncertainty in the con- 
struction of a robust LPV control law. Actuator 
faults can be parameterized as estimated fault effec- 
tiveness parameters. Thus, it is possible to formu- 
late a fault tolerance control design problem as an 
LPV control synthesis problem based on estimated 
faults parameters. 

This paper contains the following sections. In 
Section 2, an robust LPV synthesis methodology 
used in this paper is summarized. In Section 3, 
fault parameter estimation methods are presented. 
In Section 4, an LPV controller for a HiMAT ve- 
hicle is designed by using an robust LPV synthesis 
control methodology. In Section 5, the simulation 
results of the closed-loop system with an on-line es- 
timator are presented and this paper concludes with 
a brief summary in Section 6. 


2 LPV Synthesis 

2.1 Problem Statements 


where x € lZ n ‘ , € 1Z nA , e p € 1Z , y € 1Z y , 

d A e dp € n nd >, and u € n n “. All of the 

state-space matrices are of appropriate dimensions. 
An uncertainty block set A is defined as: 


A = diag{6\h , ■ • • , dnlm A n +i , • • ■ , 

:*€ fl, A t £R M ,d(A)<0 

where j3 is normalized to 1 without loss of generality. 
There exists a scaling factor set S such that 



S = {S:S> 0, SA = AS, S£Tl nAXnA }. (1) 

The input/output scaling matrices L~ l/1 and J 1/2 
are defined as 


L- 1 / 2 = 

J !/ 2 = 


S' -1 / 2 0 

0 In, 

S 1 ! 2 0 

o In.. 


g ^(nA+nd p )xfn A +nd p ) 
g ^(nA+iit p ) X (Hi + % ) _ 


The induced £ 2 -norm of a parameter dependent sta- 
ble LPV system is defined as 


l|G|| 2 «-2 


IHI-2 

vp, deC 2 , <# o IMIb 


for zero initial conditions x(0) = 0. 

The control synthesis problem is to design an 
LPV controller K{p) and scaling factor 5 in order 
to minimize the induced £2 norm of the closed-loop 
system. The robust LPV control synthesis can be 
formulated into an optimization problem: 


In this section, a control synthesis problem is de- 
fined, based on an estimated parameter vector p € 
TZ n ‘ such as actuator effectiveness. 11 Actuator effec- 
tiveness parameters represent actuator failure cases 
(actuator damage) which can be estimated using the 
estimation methods presented in Section 3. 

An LPV system can be represented as functions of 
an estimated scheduling parameter vector p with an 
uncertainty block A which captures parameter esti- 
mation errors and unmodeled dynamics. An LPV 
system can be written as: 


X 


A(p) B A (p) B p (p) B u (p) 


X 

eA 


C e (p) D AA (p) D Ap (p) D Au {p ) 


dA 

6p 


C p {p) D pA (p) D P p(p) Dp U (p) 


dp 

uV - 


Cy{p) Dy A {p) Dyp(p) 0 


_ u 


d,A — AeA, 


min || J 1 / 2 (S)Ft(P(p), K(p))L ^ (S)|| 2<-2 

* ( p), ses 

where F t (P{p), K{p)) means a lower linear fraction 
transformation (LFT). The optimization problem of 
equation (2) is not convex in K(p) and 5. When a 
estimated scheduling parameter p is fixed, the prob- 
lem is same as a D-K iteration problem. In this 
paper, the robust LPV control synthesis is formu- 
lated into two LMI optimization problems which can 
be solved using an iteration procedure. 

2.2 Control Synthesis Methodology 

In this section, an iteration procedure of solving the 
problem is presented. There is an LPV control syn- 
thesis methodology in Ref.[ 4 ] with a scaling factor. 
In Ref.[ 4 ], an LMI optimization can be formulated 
with unknown matrices X > 0, Y > 0, and scal- 
ing factor matrices J 1 < 2 and L^ 1/2 . However, an 
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LMI optimization problem in Ref.[ 4 ] has an equality 
constraint. In this paper, to avoid an equality con- 
straint, the augment LPV system with the scaling 
matrices J 1/2 and I' 1/2 is used to design an LPV 
controller. 

Suppose a scaling factor S is given. Hereafter, p 
dependence in matrices is omitted for space limita- 
tion. The augmented LPV open-loop system with 
the scaling matrices J 1 ^ 2 and L~ 1 ^ 2 can be written 


Applying the Kalman- Yakubovich-Popov (KYP) 
Lemma, 11 the LMI optimization is formulated to 
find a scaling factor 5. There exists an scaling fac- 
tor S e S which leads the induced-£ 2 norm of the 
closed-loop system is less than 7 S . The scaling fac- 
tor 5 can be determined solving the following LMI 
optimization. 

min 7 S (5) 

p>o,seS 


as: 


X 


A BiL~ 1/2 B 2 


X 

e 

= 

J l/2 Ci J 1/2 DnL~ 1/2 J 1/2 Di2 


d 

y 


Cl £>21 1/2 0 


u 


where e = [ej eJ] T and d — [ d J dJ] T . W ith as- 
sumption that £>12 and D 2 i are full column and row 
rank for all p, respectively, an LPV control synthesis 
methodology in Ref.[ 5 ’ 13 ] can be used in this paper. 

For the sake of completeness, a brief summary 
of the LPV control synthesis methodology in Ref.[ 5 ] 
is presented in this section. There exists an LPV 
controller K(p ) which leads the induced- T2 norm of 
the closed-loop system is less than constant 7. The 
LPV controller K(p ) can be constructed from the 
solution matrices, X(p) € 1Z nxn and V (p) € 1Z 
which are calculated by solving the following LMI 
optimization. 

min 7 (4) 

X(is),Y {p) 


subject to 


A t P + PA + P PBi cl PB 2cl 
Bj cl P S 0 

Bl,P 0 -I _ 


V 

O 

hi 

O 
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SCi cl SDu ci SDu c , 

4- 

dTi cI m i„ 


_7 _l C , 2 c( Y~ lD 22 c ,_ 
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The iteration procedure to solve the problem in 
equation (2) is follows: 


1. Design an LPV controller K (p) for a system 
from the LMI optimization in equation (4) with 
fixed 5. At the first iteration, S is assumed as 

1. 

2. Solve the LMI optimization problem in equa- 
tion (5) over P(p) and S based on the closed- 
loop system with the designed LPV controller 


subject to the LMI constraints defined in Ref.[ 5 ]. 
The realization of an LPV controller from the so- 
lution matrices X and Y are taken from Ref.[ 5 ]. 

The benefit of the LPV synthesis methodology 
is that there is no limitation of an affine functional 
form of LPV system state-space matrices. Since the 
LMI constraints are evaluated at grid points over all 
scheduling parameter spaces, an LPV system should 
be just a continuous and differentiable function of a 
scheduling parameter. The disadvantage of the LPV 
control synthesis methodology is that robust stabil- 
ity over all parameter spaces is not guaranteed unless 
choosing appropriate number of grid points. 

Suppose there exists an designed LPV controller 
K(p) which stabilizes the augmented LPV system. 
The closed-loop LPV system with a given controller 
is written as: 


Xcl 


A cl 


B 2 cl 



= 

Ci cl 

Due 

B>2 l c . 




Ci cl 

Dn cl 

D 22 cI \ 

L 


where x c i = [ x T x -[] T , B\ cl 

C Vcl = S 1/2 C Ul , D Ucl = S'' 2 Dn cl 
S l / 2 D-n cl , and D 2 le( = D 2lcl S~^ 2 . 


, d /x — j 


S- 1 ' 2 , d 21cI = 


K(p). 


3. Generate an augmented LPV system with the 
scaling factor S 


G t +i(p) 


gl/2 

0 



S-l/2 

0 


0 

I 


( 6 ) 


where G t (p ) is an LPV model at the i th itera- 
tion. 


4. Iterate over step 1 to 3 until convergence or ter- 
minate iteration based on satisfaction with a 
designed LPV controller. 


The iteration method can not guarantee finding 
global solutions of K and S since the problem in 
equation (2) is not convex in K and 5. Also, there is 
no guarantee of convergence in the iteration process. 
When a designed controller is fixed to calculate the 
scaling factor 5, the matrix P{p) can be calculated 
from the solution matrices X and Y and fixed in 
the LMI optimization. 5 Note that the matrix P(p) 
is also set as an unknown matrix in the LMI opti- 
mization in equation (5) in this paper, to relax the 
constraints of fixing the LPV controller K (p). 
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Bias-free state estimator 


3 Parameter Estimation 

This section briefly describes the formulation of a 
real parameter estimation problem, which, when 
specialized to the actuator effectiveness estimation, 
transforms a fault (loss of actuator effectiveness) tol- 
erant control problem to a robust LPV control prob- 
lem. The development of this section follows that in 
Ref.[ u ]. 

The estimator is based on a linear discrete design 
model of the form: 

' ' 

Ifc+l = Mi x k + [&l7fc ‘ ' ‘ ] 


+B%Uk + U>k 

= Aix k + Efak+Biu k + w x k , (8) 

7 fc+i = Ik Awl, ( 9 ) 

y k = C d k x k +v k , ( 10 ) 


where x k £ R n x , 7 k £ R n “ 7 u k € R Uu and y k £ R v 
are the state, bias, input, and output variables, re- 
spectively. The discrete model can be obtained from 
a continuous model via, for example, the Euler s rule 
with a sampling period T s , which preserves the func- 
tional dependence of the u B ” matrix on 7. The 
bias vector 7 with component -1 < 7* £ 0 relates 
to a actuator failure parameter. It is obvious that 
E% = Bfxdiag{u l k ,--‘ ,u£ u }- ^ w k & ndv k denote 

the white noise sequences of uncorrelated Gaussian 
random vectors with zero means and covariance ma- 
trices Q%, Ql and R k , respectively. 

The minimum variance solution is obtained by a 
direct application of the two-stage Kalman filter al- 
gorithm of Keller and Darouach, 1 ^ with constant co- 
efficient matrices in Ref.[ 14 ] replaced by time- varying 
matrices. The filter is decoupled into four sets of 
equations. They given as follows. 

— Optimal bias estimator 


7k+l\k 

-sc 

-sc 

c- 

II 

P 1 

r k+l\k 

= p k\ k + Ql 

7*+i|jfe+i 

— 7fc+i|* + ^b+i^fc+i 

K 2+i 

_ p7 fiT 

~ r k+i\k n k+i\k 


(Hk+\\kPk+i\kH k +l\k d - Sk+l) 5 

P 1 

Gt+ijfc+i 

— {I ~ Kk+lHk+l\k)Pk+l\k‘ 


%k+l\k 

= A k x k \ k + B k u k + W k l k \k ~ U+i|fc7fc|* 

px 

r fc+ljJt 

= AlPfa ( A d f + Ql + W k P^ k Wj 


~Vk+l\kP k+ i\ k ^ k+l\k 1 

X k +i\k+\ 

= £*;+l|fc + K k+ \{Hk + l ~ Cfc+l^*+l|*)’ 

Kk+ 1 

= P k+ \\ k {Ck+\) T {Ct+iPk+i^Ct+i) 


+Rk+l}'\ 

^fe+l|fc + l 

II 

1 

** H 

§ 
+ N 


where the filter residual and its covariance are given 
as 


rk+i — Vk+i ~ 

Sk + 1 — + Rfe+l. 

— Coupling equations 

W k = A d k V k]k +E d k , 

u + i, fc = w k p^ k (pz +l lk r l , 

Hk+i\k = C k+ ]V k+ i\ k , 

U-t-L|fc+l = U+l|* — R k + l^k+l\k- 

— And finally the compensated state and error co- 
variance estimates 

£fc+i|fc+i = ®*+i|Jfe+i + U+i|fc+i7*+i|fc+i> 

^)t+i|fc+i = -Pjfc+i|fe+i + U+i|fc+i^fc+i|fc+iU+i|fc-H‘ 

A further measure is taken to modify the above fil- 
tering algorithm so that the estimates become more 
responsive to abrupt changes in the control effective- 
ness factors. 

A well known technique for estimating time- 
varying parameters is the use of forgetting factors. 
The basic idea is to enable a recursive algorithm to 
discount the past information so that the filter is 
more apt to recognize the changes in the system. 
Since the time update of the bias estimate governed 
t>y 7fc.j_i|A: — 7fc|fc is the dominant opposing force to 
acknowledge the abrupt changes in the biases, for- 
getting factors introduced into the time propagation 
equation P^ l]k = P^ k + Ql of the bias covariance 

is likely to function most effectively. 

Assume that covariance P^ k adequately describes 
the bias estimation error along both temporal and 
spacial directions under the normal system oper- 
ation condition. Then this covariance provides a 
basis for the selection of forgetting factors. The 
bias estimates should be prevented from being im- 
petuous, as well as from being indifferent to the 
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changes shown in the measurements. A technique 
suggested in Ref. [^] amounts to select forgetting 
factors that would force the adjusted covariance in 
p 1 , = P ?, . + 07 to stay within some prescribed 

bounds 

& min ^ 1 P^y\ k - (H) 

where <j m i n and a ma x are positive constants with 0 < 
(7 min < cr max < oo, and I is the identity matrix. Let 
the dyadic expansion of P^ k be given by 


Ai* = X>L e *(4A < 12 ) 

i- 1 

where a \ [k , .... a^ k are the eigenvalues of P^ k with 
c^ |fc > ... > a”“ fc , and e 1 ,...^ are the correspond- 
ing eigenvectors with ||ej. || = ... = 1 1 ^ ^ u || = 1- Equa- 
tion (12) can then be expressed as 


P 1 


= E : 


i— 1 


#el(4) T + QZ, 

*L- 


0< A' < 1. 


Following the argument in Ref.[ lr5 ], the forgetting 
factor A^ can be chosen as a decreasing function of 
the amount of information received in the direction 
ej.. Since eigenvalue aj.^ of P^ k is a measure of 

the uncertainty in the direction of e l k , a choice of 
forgetting factor Aj^ based on the above constraints 
can be 


canards two outputs: angle of attack a in radians 
and pitch angle 0 in radians; and four states: velocity 
V in ft /sec, angle of attack, pitch rate q in rad/ sec, 
and pitch angle. The open-loop model is 


r 


r 


* “ 


- 
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A B 


X 


<5e 






1 ^ 


y. 


1 

O 


u 


A 


where 



-0.0226 

-36.6 

-18.9 

-32.1 



0 

-1.9 

0.98 

0 



0.012 

-11.7 

-2.63 

0 



: 0 

0 


1 


0 



0 

0 







-0.414 

0 


— 

0 

1 0 

0 


-77.8 

22.4 



0 

0 0 

1 


0 

0 







A failure parameter vector r = [n 72 ] t is deter- 
mined by the two actuator effectiveness parameters 
n and t <2 of elevon and canard actuators, respec- 
tively. Assume that the failure parameters linearly 
enter in the model. The state-space model of the 
HiMAT vehicle is written as 


Al = 


Ao , 

a \\k [ a min + 




Q A:|fc > a niax, 
a lc|A: < Omax- 

(13) 


The estimation algorithm discussed in this section 
will be seen to have been applied successfully to a 
HiMAT model. 


4 Control of HiMAT Vehicle 


x — Ax + B(t)u, y = Cx, 

where .4 and C are constant matrices and B(t) 
[b\Ti l> 2 t~ 2 } ■ The vectors b\ and bo are the columns 
of B. The actuator failure parameters can be es- 
timated using the estimation method described in 
Section 3. However, there is estimation error S T — 
[j T| 5 T2 \ t . The actuator failure parameter vector r 

is written as 

r = f + 6 t , 


In this section, the robust LPV control synthesis 
methodology described in Section 2 is applied to 
control a HiMAT vehicle for actuator failure cases. 
Recall it is assumed that actuators are failed one at 
a time. Thus, the control reconfigurability of the 
HiMAT vehicle never goes to zero over the entire 
failure cases. 16 The system variations due to actu- 
ator failures can be modeled as an LP\ system, a 
function of an estimated scheduling parameter. 

4.1 LPV Model of HiMAT Vehicle 

The model of the HiMAT vehicle taken from the /i- 
synthesis Toolbox 12 has two inputs: elevons <5 e and 


where f is an estimated value. The estimated actu- 
ator failure parameter can be written as: 


f, =l+7 1 , 

f-2 = 1 + 7 2 - 


where 7 1 and 7 2 are bias parameters in equat ion (9 ). 
The estimation error bound is assumed as \fo^6 T < 
0.05 for each actuator failure case. The matrix B(t) 
is rewritten as: 


B(t) = B 


fi + 0.05.5! 
0 


0 

T2 4- 0.05<52 


(15) 
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where the real uncertainty parameters <5i and d- 2 vary 
from -1 to 1, respectively. 

The LPV model of the HiMAT vehicle is 


x = Ax + B 


T\ 0 

0 T‘2 


u + B 


0.05 0 

0 0.05 


z — u, y — Cx , 


w, (16) 
(17) 


w = A z, 



0 


(18) 


The LPV model in equation (16) is a linear function 
of a parameter vector f. However, the parameter 
vector t can not be chosen as a scheduling param 
eter since we consider that the actuators are failed 
one at a time. Thus, f\ and f2 can not be zero simul- 
taneously. To describe the failure cases, a synthetic 
scheduling parameter p is introduced as 0 < p < 2. 


0 < fi < 1, f 2 = 1 

fi = 1, f 2 = 1 

T\ — 1 , 0 < f‘2 < 1 

Note that the LPV model of the HiMAT vehicle is 
not an affine function of a scheduling parameter p . 


0 < p < 1 

p = 1 

1 <p<2 


4.2 LPV Controller Design 

The control objective is to track a pitch angle com- 
mand at actuator failure cases. A designed LPV 
controller should robustly stabilize the HiMAT ve- 
hicle over the failure parameter variations. The con- 
troller synthesis problem is formulated as a model 
matching problem in Figure 1. The ideal response 



Figure 1: Interconnection structure for the model 
matching problem. 


//-synthesis Toolbox. 12 The performance weighting 
function W p and unmodeled dynamics \\ n axe also 
taken from the example in //-synthesis Toolbox. 12 
The sensor noise is modeled as white noise with 0.6 
amplitude for angle of attack and pitch angle mea- 
surements. The weighting functions in Figure 1 are 


Ti 


= 40 


1 

s/0.8 + 1 1 
s /50+1 


= 0.2 


s/ 0.05 + r 
s /5 + l J 


s/1000+1 


2x2i 


w noise 


= 0.01/2x2- 


The control synthesis problem of the HiMAT vehi- 
cle is formulated to minimize the induced- C 2 norm 
of the augmented LPV system with the weighting 
functions. 

To solve the control synthesis LMI optimization 
problem in equation (4), basis functions for X (p) 
and Y{p) are required since X and Y are assumed 
as functions of p. 


X (p) 

* (20) 
Y(P) =Yj9iWi, YjZK n * n . 
j 

where basis functions fi{p) and gj{p) are given be- 
fore solving the LMI optimization in equation (4) 
over Xi and Y v There is no analytic method to 
choose optimal basis functions for X and V in gen- 
eral. In this paper, the basis function set is defined 
as {1,1 Ip, p) for X and Y to help the LMI optimiza- 
tion for total failure cases (fi = 0,f2 =0), since the 
functions }i{p) and gj(p) are related with sensitivity 
of unknown matrices Xi and Yj, respectively. Note 
that it is not necessary to define that gj(p) is equal 

to ft(p). _ 

Since X and Y are functions of p, the parameter 
rate bound is required to solve the LMI optimization 
in equation (4). Recall that the scheduling param- 
eter is an actuator failure parameter. Thus, for ex- 
ample, the scheduling parameter can suddenly vary 
from 1 (no failure case) to 2 (total canard failure 
case). In this paper, the parameter rate bound is 
assumed as |^| < 100 to capture sudden variations 
of the scheduling parameter. 

To make the LMI optimization computationally 
tractable, the LMI constraints are evaluated at the 
following grid points: 


model Ti of pitch angle is taken form the example in 


pe {p|0.01, 0.1, 0.2, ••• , 1.9, 2}. 
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Iteration 

7 

S 

1 

1.23 

diag([l, 1,1,1]) 

2 

0.71 

diag( [0.497, 0.168, 1.186, 1.277]) 

3 

0.60 

diag( [1.404, 1.289, 1.801, 1.453]) 

4 

0.54 

diag( [1.580, 2.193, 2.100, 1.187]) 

5 

0.85 

diag([2.007, 1.430, 2.506, 1.849]) 


Table 1: 7 values in the LMI optimization in equa- 
tion (4). 

Also, the same grid points are used to solve the LMI 
optimization in equation (5). Since the matrix P is 
related with X and Y\ the basis function set for P 
is chosen as {1, 1/p, p}. 

In this paper, the scaling factor 5 is assumed as 
constant over all scheduling parameter variations. 
The 7 values and the scaling factor 5 for each it- 
eration are written in Table 1. The scaling factor S 
is associated with the uncertainty block A which is 

A = diag([S 1, 5 $ele, <5can])* 

Recall that the real uncertainty parameters <$1 and 
82 are associated with elevon and canard actuator 
failure parameters, respectively. The multiplicative 
uncertainty parameters S e i € and S can are also associ- 
ated with elevon and canard control channels in Fig- 
ure 1. The iteration process is stopped at the b th 
iteration since the 7 value at the iteration is greater 
than the previous iteration. Recall that the iteration 
process is not guaranteed to be converged. However, 
the performance index 7 in the LMI optimization of 
equation (4) is significantly reduced from 1.23 to 0.54 
by using the scaling factor S. In the remainder of 
this paper, the LPV controller for the HiMAT vehi- 
cle denotes the designed LPV controller at the 4 th 
iteration. 


5 Simulations 

In this section, the designed LPV controller is ap- 
plied to control the HiMAT vehicle for actuator fail- 
ure cases. The HiMAT vehicle is simulated at a pre- 
defined failure scenario. One of the failure scenarios 
is defined as: the canards are failed from 1 to 10 


seconds and the elevons are failed from 20 to 40 sec- 
onds. The actuator failure parameter for the failure 
scenario is set as: 

{ 0.01, 1 <t< lOsec, 

1.00, 0 <t< Isec, 10<£< 20sec 

2.00, 20 < t < 40sec. 

(21) 

For comparison, the LPV controller for the HiMAT 
vehicle is simulated with and without actuator fail- 
ures. The simulation results are shown in Figure 2. 
The pitch angle commands are given as 10° at 1 sec, 
0° at 10 sec, and 10° at 20 sec, sequentially. In this 



time (sec) 

Figure 2: The LPV controller simulations with and 
without actuator failures. 

simulation, we assumed that the failure parameter 
is measured as true value shown in the bottom plot 
in Figure 2. It is observed that the LPV controller 
achieves the desired performance level of tracking 
pitch commands at actuator failure cases. From the 
control signals for actuators in Figure 2, the LPV 
controller significantly uses the canard actuator to 
track the pitch commands when the elevon actuator 
is failed from 20 sec to 40 sec. Also, the LPV con- 
troller keeps the elevon actuator signals close to zero 
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at the elevon actuator failure case. 

For the robustness of the LPV controller, the un- 
certainty parameters &\ and £2 are chosen as certain 
numbers such that = diag([8\, <5 2 ]) T ||A P || < 1 

and are folded into the LPV model of the HiMAT 
vehicle. The LPV controller is simulated with the 
perturbed system at the same failure scenario de- 
scribed in equation (21). Two examples of uncer- 
tainty parameters are shown in Figure 3 for space 
limitation, “pertl” and “pert2” in Figure 3 denote 
[Si S 2 ] = [1 1] and [Si S 2 ] = [-1 - 1]» respec- 

tively. The simulation results show that the LPV 
controller can robustly stabilize the perturbed sys- 
tem and achieve the desired performance level of 
tracking the pitch angle commands. 

For comparison, an H 0 0 controller as a linear time 
invariant controller for the HiMAT vehicle is de- 
signed at p = 1 (without failure cases) with the same 
weighting functions described in Section 4, using fi- 
synthesis Toolbox. 12 The H controller is applied 
for the actuator failure scenario of equation (21). 
The bottom plots of Figure 3 show that the H 0 0 con- 





Figure 3: Time responses of pitch angle with LPV 
and Hoc controllers. 

troller can achieve the desired performance level at a 
canard failure case(l < t < 10 sec) even though the 


controller is designed at p = 1. It shows that the 
Hoc controller inherently has some level of robust- 
ness. However, the controller can’t achieve the 
desired performance level at the elevon failure case. 
Since the canards are less effective in controlling the 
pitch movement in comparison with elevons, 16 the 
elevon failure can significantly affect on tracking the 
pitch commands. 

The control signals of the LPV and #00 con- 
trollers are shown in Figure 4 at the same failure 
scenario. Since the control signals of the LPV con- 




Figure 4: Control signals for each actuator with the 
LPV and H 0 c controllers. 

troller with/ without perturbations of S 1 and 62 are 
same, the “pertl” and “pert2” simulation results are 
plotted in Figure 4. It is observable from Figure 4 
that the elevon signals of the H 0 0 controller are sig- 
nificantly increased at the elevon failure case. That 
is not appropriate to implement the Hoo controller 
into the real system. However, the LPV controller 
significantly reduces the elevon signals and increases 
the canard signals to compensate the elevon failure 
in Figure 4. 

Now, the fault parameter is estimated with the 
on-line estimator integrated with the LPV controller 
shown in Figure 5. The on-line estimator in Figure 5 
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has two parts: one is the two-stage discrete Kalman 
filter and the other is a simple logic in equation (19) 
to calculate the fault parameter from the bias esti- 
mates. To simulate the Kalman filter with the Hi- 
MAT vehicle, several parameters should be defined. 
In this paper, sampling time is set as 0.01 sec to 



Figure 5: Simulation block diagram of a closed-loop 
system with an on-line estimator. 

capture the open-loop dynamics of the vehicle. The 
covariance matrices Q x k , Ql and Rl described in in 
Section 3 are set as constant matrices: 

Ql — 3diag([l, 0.01 2 , 0.01 2 , 0.01 2 ]) 

Ql = 3dm#([0.05 2 , 0.05 2 ]) 

Rl = 3dias([0.01 2 , 0.01 2 ]) 

since the covariance matrices are related with con- 
vergence of the estimator and noise allowance con- 
sidered in the control synthesis in Section 4. The 
initial values of estimated state x 0 |o and biases 7 0 | 0 
are set as [0 0 0 0] T and [0 0] r with assumption of 
perfect estimation at t = 0. The initial covariance 
matrices P^ 0 and P^ 0 are defined as 10/ 2 and 10/ 4 . 
We choose the different values such as 0, 100 for the 
initial covariance matrices for simulations. The bias 
estimates are not dependent on the initial values. 

The most delicate parts of the bias estimation 
mechanism used in this paper are to define values of 
Ao, Qrmn and a m ax in equation (13) for covariance 
propagation. There is no a analytic way to define 
the values in general. Also, it is not known how the 
values affect on robustness of the closed-loop system. 
In this paper, we simulate the on-line bias estimator 


Case 

Ao 

Ckmin 

®max 

1 

0.90 

10 

100 

2 

0.95 

10 

100 

3 

1.00 

10 

100 

4 

1.00 

10 6 

10 7 


Table 2: The different sets of the varying forgetting 
factors 


based on the different sets of the values in Table 2. 
The fault parameter on-line estimate results are 




Figure 6: Fault parameter estimates according to 
the cases in Table 2. 

shown in Figure 6. It is noticeable that bias esti- 
mates are dependent on the cases. From the top 
plots of Figure 6, it is noticed that the value of A 0 
affect on convergence rate of estimation. When A 0 
is set as 1 at Case 3, the bias estimate is not con- 
verged at canard failure situation. For Case 3, the 
fault parameter estimator returns false failure indi- 
cation which is elevon failure instead of canard fail- 
ure. Since the canards of the HiMAT vehicle have 
less controllability than the elevons do, it is much 
harder to detect the canard failure than to do the 
elevon failure. In Case 1,2, and 4, the false failure 
indication around 1 sec is shown in Figure 6. For 
Case 4, the covariance matrix P y is immediately 
high value after one step integration since a mrn is 
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elevons(dey) <» (deg) 


defined as 10 6 . Case 4 is not realistic. However, it 
is easily noticed how P 1 affects on the estimator. 

Simulation results at Case 1 are shown in Fig- 
ure 7. “TV” and “VFF’ in Figure 7 denote that 
the LPV controller is evaluated at the true values 
(TV) of the failure parameter and at the estimated 
failure parameter with a variable forgetting factor 
(VFF), respectively. It is noticeable that the LPV 
controller evaluated at the estimated parameter can 
achieve the desired performance of tracking the pitch 
commands. The difference between estimated failure 
parameter and true failure parameter is very small at 
a steady state. However, there is a delay to estimate 
the failure parameter, which is not included in the 
LPV control synthesis process. The delay is help- 




tme{sec) 


Figure 7: Simulations with the on-line estimator. 


6 Conclusion 


In this paper, the LPV controller is designed based 
on the estimated scheduling parameter. The system 
variations due to actuator failures are modeled as 
functions of estimated parameters and bounded pa- 
rameter estimation errors represented by LFT form 
with the uncertainty block. An LPV controller syn- 
thesis problem with bounded parameter errors is for- 
mulated into the two LMI optimizations which are 
solved by the iteration method. The performance 
level of the closed-loop system with the designed 
LPV controller is reduced by the constant scaling 
factor on the uncertainty block. 

The iteration approach of the LPV synthesis 
methodology is applied to control of the HiMAT ve- 
hicle at actuator failure cases. It is assumed that the 
actuators are failed one at a time. Actuator failure 
parameters of the HIMAT vehicle are estimated as 
biases using an augmented Kalman filter. The LP\ 
controller evaluated at the estimated failure parame- 
ter is simulated with the HIMAT vehicle which varies 
as true values of the failure parameter. The simu- 
lation results show that the LPV controller achieves 
the desired performance level of tracking pitch angle 
commands for actuator failure cases. 
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ful to make the scheduling parameter rate bounds 
which is one of the assumptions of the LPV control 
synthesis methodology. However, the large delay of 
fault parameter estimations is not appropriated for 
a fault tolerant control system. It is not studied 
here to include the delay term in the LPV control 
synthesis methodology. It is unknown to find out 
the optimal delay for reconfigurability of the LPV 
controller in the fault tolerant control system. 


References 

[1] Blanke, M., Izadi-Zamanabadi, R., Bogin 
S.A., and Lunau, C.P., “Fault- Tolerant Control 
System- A Holistic View,” Control Engineering 
Control Practice , Vol. 5, No. 5, 1997, pp. 693- 
702. 


10 

American Institute of Aeronautics and Astronautics 






[2] Chen, J., Patton, R.J., and Chen, Z, “An 
LMI approach to Fault-Tolerant Control of 
Uncertain System,” in Proceeding of the 
1999 IEEE/ ISIC/CIR A/ISAS Joint Confer- 
ence, 1998, pp. 175-180. 

[3] Ganguli, S., Marcos, A., and Balas, G., “Recon- 
figurable LPV Control Design for Boeing 747- 
100/200 Longitudinal Axis,” in Proceedings of 
American Control Conference, 2002, pp. 3612 
3617. 

[4] Apkarian, P. and Adams, R., “Advanced Gain- 
Scheduling Techniques for Uncertain Systems,” 
IEEE Transactions on Control Systems Tech- 
nology, Vol. 6, No. 1, 1998, pp. 21-32. 

[5] Wu, F., Control of Linear Parameter Vary- 
ing Systems. PhD thesis, Department of Me- 
chanical Engineering, University of California, 
Berkeley, 1995. 

[6] Wu, F., Packard, A., and Balas, G., “LPV 
Control Design for Pitch- Axis Missile Autopi- 
lots,” in Proceedings of the 34th IEEE Confer- 
ence on Decision and Control, (New Orleans, 
LA), 1995, pp. 188-193. 

[7] Shamma, J., “Gain-Scheduled Missile Autopi- 
lot Design Using Linear Parameter Varying 
Transformations,” Journal of Guidance, Con- 
trol, and Dynamics, Vol. 16, No. 2, 1993, 

pp. 256 - 261 . 

[8] Balas, G., Fialho, I., Packard, A., Renfrew, J., 
and Mullaney, C., “On the Design of LPV Con- 
trollers for the F-14 Aircraft Lateral-Directional 
Axis During Powered Approach,” in Proceed- 
ings of the American Control Conference, (Al- 
buquerque, NM), 1997, pp. 123-127. 

[9] Balas, G., Ryan, J., Shin, J-Y., and Garrard, 
W., “A New Technique for Design of Con- 
trollers for Turbofan Engines,” in AIAA 34th 
Join Propulsion Conference, (Cleveland, OH), 
1998, pp. 1-6. AIAA-98-3751. 

[10] Shin, J.-Y., Worst-Case Analysis and Linear 
Parameter- Varying Gain-Scheduled Control of 
Aerospace Systems. PhD thesis, Department of 
Aerospace Engineering and Mechanics, Univer- 
sity of Minnesota, 2000. 

[11] Wu, N. Eva, Zhang, Y., and Zhou, K., “Detec- 
tion, Estimation, and Accommodation of Loss 
of Control Effectiveness,” International Jour- 
nal of Adaptive Control and Signal Processing, 
Vol. 14, 2000, pp. 775-795. 


[12] Balas, G., Doyle, J., Glover, K„ and Packard, 
A., p Analysis and Synthesis Toolbox. Natick, 
MA: The Mathworks, Inc., 1995. 

[13] Becker, G., Quadratic Stability and Perfor- 
mance of Linear Parameter Dependent Sys- 
tems. PhD thesis, Department of Mechanical 
Engineering, University of California, Berkeley, 
1993. 

[14] Keller, J. K. and Darouach, M., “Optimal two- 
stage Kalman filter in the presence of random 
bias,” Automation, Vol. 33, 1997, pp. 1745- 
1748. 

[15] Parkum, J. E., Poulsen, N. K., and Holst, 
J., “Recursive forgetting algorithms,” Inter- 
national Journal of Control, Col. 55, 1992, 
pp. 109-128. 

[16] Wu, N. Eva, Zhou, K., and Saloman, G., “Con- 
trol Reconfigurability of Linear Time-Invariant 
Systems,” Automatica, Vol. 36, 2000, pp. 1767- 
1771. 


11 

American Institute of Aeronautics and Astronautics 




AIAA Guidance, Navigation, and Control Conference and Exhibit 
5-8 August 2002, Monterey, California 


AIAA 2002-4546 


LINEAR PARAMETER VARYING CONTROL SYNTHESIS 
FOR ACTUATOR FAILURE, BASED ON ESTIMATED 

PARAMETER * 
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Abstract 

A robust linear parameter varying (LPV) control 
synthesis is carried out for an HiMAT vehicle subject 
to loss of control effectiveness. The scheduling pa- 
rameter is selected to be a function of the estimates 
of the control effectiveness factors. The estimates 
are provided on-line by a two-stage Kalinan estima- 
tor. The inherent conservatism of the LPV design 
is reduced through the use of a scaling factor on the 
uncertainty block that represents the estimation er- 
rors of the effectiveness factors. Simulations of the 
controlled system with the on-line estimator show 
that a superior fault- tolerance can be achieved. 


1. Introduction 

One of the most effective methods for achieving sin- 
gle aircraft accident prevention 1 is to design flight 
control laws that are fault-tolerant. An active fault- 
tolerant control (FTC) system requires its control 
law to react to actuator/ sensor faults through re- 
configuration. It is important that the design of a 
fault detection and isolation (FDI) mechanism in an 
FTC system be properly integrated with the design 
of a fault tolerant control mechanism. 2 

Fault-tolerant control designs for aerospace vehi- 
cles based on linear matrix inequality (LMI) opti- 
mization solutions can be found in several reports. 
For example, in Ref.[ 3,4 ], fault tolerant controllers 

* Copyright © 2002 by J-Y. Shin, N. E. Wu, and C. Bel- 
castro. Published by the American Insitute of Aeronautics 
and Astronautics, Inc. with permission. 

t Staff Scientist, ICASE, NASA Langley Research Center, 
Hampton, VA 23681. AIAA Member 

t Professor, Department of Electrical Engineering, Bing- 
hamton State University of New York, Binghamton, New 
York, 13902. 

S Researcher, Guidance Control Brench, NASA Langley 
Research Center, Hampton, VA 23681. 


are designed based on information provided by sep- 
arate FDI modules that are assumed to have iden- 
tified faults, while faults are modeled as functions 
of model parameters. In particular, when faults 
are modeled as a set of varying parameters in a 
linear parameter varying (LPV) system, a known 
LPV control synthesis method 5 ’ 6 can be attempted. 
The synthesis is formulated into an LMI optimiza- 
tion problem based on the LPV system whose state- 
space matrices are functions of a scheduling param- 
eter vector. A scheduling parameter vector must be 
measurable in real time and its rate of change is re- 
quired to be bounded. 5,6 It is possible that an LPV 
controller can robustly stabilize a close-loop system 
and achieve a desired performance level over the en- 
tire parameter space. The LPV synthesis methodol- 
ogy has been successfully applied to synthesize con- 
trollers for the pitch-axis missile autopilots, 7,8 the F- 
14 aircraft lateral-directional axis using powered ap- 
proach, 9 turbofan engines 10 and the F-16 aircraft. 11 

In this paper, fault parameters are considered 
as scheduling parameters in LPV control synthesis, 
whose rate of change is defined as a sufficiently large 
number but not infinite. In practical, the bounded 
large number of scheduling parameter rate of change 
can allow abrupt fault occurrence in a reconfigurable 
LPV controller in Ref.[ 4 ]. Since the fault parameters 
cannot be directly measured via sensors in real time, 
they are estimated on-line instead using a two-stage 
adaptive Kalman filtering scheme. 12 To deal with 
the uncertainty due to estimation errors which al- 
ways exist in practice, in this paper, modeling of 
fault parameter estimation error is also included as 
the part of the problem formulation in the form of 
an uncertainty block, and this estimation error en- 
ters in the construction of an LPV controller. Thus, 
the LPV controller in an FTC system can allow to be 
easily integrated with an FDI system-the on-line pa- 
rameter estimator. The LPV control synthesis prob- 
lem can be formulated into a robust LPV control 
synthesis. 5 In this paper, the robust LPV control 
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